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A model for GUT baryo genesis, coherent baryogenesis within the framework of super- 
symmetric SO(10), is considered. In particular, we discuss the Barr-Raby model, where at 
_ ' the end of hybrid inflation charge asymmetries can be created through the time-dependent 

higgsino-gaugino mixing mass matrix. These asymmetries are processed to Standard Model 
^ . matter through decays via nonrenormalizable (B — i)-violating operators. We find that a 

baryon asymmetry in accordance with observation can be generated. An appendix is devoted 
t— I ' to provide useful formulas and concrete examples for calculations within SO(10). 
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Grand Unified Theories (GUTs) generically predict a scalar potential and thereby a large 

Oh! 

amount of vacuum energy when the scalar fields are displaced from the minimum which breaks 
the symmetry down to the Standard Model. It is hence often argued that cosmic inflation may be 

• • • ' n 

J> , implemented by the scalar field dynamics of a GUT [1|. A well-known paradigm is super symmetric 

X ■ (SUSY) hybrid inflation , and it has led to various models using different grand unified gauge 

groups, see e.g. [4, 0,0]. 

Another feature of GUTs, which is of possible relevance for cosmology, is obviously the violation 
of baryon minus lepton number, B — L, due to the unification and mixing of baryons and leptons, 
since this can lead to mechanisms for generating the observed baryon asymmetry of the Universe 
(BAU). While leptogenesis is strictly speaking not necessarily implemented into a GUT - it can be 
operative within the Standard Model minimally extended by right handed neutrinos with Majorana 
mass terms - we suggested a scenario relying on baryondepton unification, coherent baryogenesis 2I 
and implemented it within a Pati-Salam supersymmetric hybrid inflationary model. We give 
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a brief review of this mechanism and our calculational formalism inspired by kinetic theory in 
section El 

Being a product group, the Pati-Salam gauge group 9] gives, strictly speaking, not rise to a 
GUT, and it is therefore desirable to devise models based on the proper GUT SO(10). A way 
of breaking the SO(10)-symmetry, which is particularly suitable for SUSY hybrid inflation, is the 
mechanism proposed by Barr and Raby. Kyae and Shah extended this model by global symmetries 
which restrict the superpotential to contain only couplings consistent with hybrid inflation HJ. 
This model has however a rather complicated Higgs sector. 

The purpose of this paper is to show that coherent baryogenesis naturally occurs within a 
SUSY-SO(10)-framework. Therefore, we want to keep the discussion as simple as possible and use 
the minimal superpotential suggested by Barr and Raby, leaving aside the issues which Kyae and 
Shah focus on. In turn, we also point out that we do not extend the minimal model by ad hoc 
terms just in order to make our mechanism viable. 

In section |3] we put together the features of the Barr- Raby model which are important for our 
baryogenesis scenario and present in some detail the derivation of the higgsino-gaugino mixing mass 
matrix. We also devote an appendix to the conventions and techniques we apply for calculations 
within SO(10), with the intention to make this article self-contained and easily comprehensible to 
the reader who is not familiar with SO(10)-model building, and furthermore we want to provide a 
useful help for accessing the papers mentioned above. 

Putting our considerations into work, we present a numerical study of SO(10)-coherent baryoge- 
nesis in section|lJ There, we also discuss the decay processes of the higgsinos via nonrenormalizable, 
{B — L)-violating couplings. The result is an estimate of the produced baryon asymmetry for a 
particular set of parameters. 



2. COHERENT BARYOGENESIS 



Coherent baryogenesis relies on the production of particles due to a time dependent mass 



term, a phenomenon which we refer to as 



for scalars 



llfl, 11, 12, 3 



and for fermions 



preheating. Preheating has been extensively studied 



using the technique of Bogolyubov transfor- 



mations. A special feature of fermionic preheating is that fermions can be amply produced when 
their mass term crosses zero. When one considers instead of a mass term a mixing mass matrix, 
charge C and parity P may be violated and an asymmetry can be stored within different fermionic 
species. A nonsymmetric mass matrix however leads to the violation of the orthogonality of par- 
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ticle and antiparticle modes. One therefore needs a formalism which is independent of a basis in 
terms of particle and antiparticle creation and annihilation operators and thereby generalizes the 
Bogolyubov transformation approach. This is developed in Ref. IT] and shall be briefly reviewed 
in the following. 

We consider several fermionic flavours, mixing through a mass matrix M{rj), which is a func- 
tion of the conformal time rj. In a spatially flat Friedmann-Lemaitre-Robertson- Walker Universe, 
described by the metric g^ v = a 2 {r]) x diag(l, —1, —1, —1), we rescale the fields such that for the 
mass terms, there is the replacement M(rj) — * a{rj)M{rj). 

Our goal is to compute the charge density, which is a bilinear form in the fermionic fields. We 
therefore introduce the Wigner function 

iS<j(k,x) = -J dVe ifc ' r (^(x-r/2)^(x + r/2)>, (1) 

where i, j are flavour indices and (i7°S' < )^ = i7°S' < is hermitean. The Wigner transform is the 
Fourier transformation of the two point function w.r.t. its relative coordinate while keeping the 
center of mass coordinate fixed. One can hence consider it as an analogue to a classical phase space 
density defined in quantum theory. Since we assume here spatial homogeneous conditions, one can 
ignore the center of mass coordinate in the following and consider iS < as a Fourier transform. 
Our formalism is applicable for a 2-point function with general density matrix, but in view of our 
applications in inflation we prefer to write it with respect to the vacuum from the outset. 
When decomposing the mass matrix Af into its hermitean and antihermitean parts, 

M H = 1 (M + Aft) , Af4 = | (M - Aft) , 
we find that \S < obeys the Wigner space Dirac equation 

(t + ~ ( M h + h 5 M A )e-& 9 «°\iS< = . (2) 

The mass matrix M emerges generically from Yukawa couplings to scalar field condensates, 
Cyu = —y<ffl>RipL + h- c - I n t ne model we consider here, Af is the higgsino-gaugino mixing mass 
matrix. 

A crucial point is the time-dependence of Af , which is not only the source of particle production. 
The matrices Af and dM/dr] both contribute to CP-violating phases, which - provided Af and 
dM/dr) are linearly independent - can not be removed by time-independent redefinitions of the 
fermionic fields. 
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In order to simplify the Wigner-Dirac equation P|). which is, besides the flavour indices, also 
endowed with a 4 x 4 spinor structure, we make use of the fact that for spatially homogeneous 
ijoS^, the helicity operator h = k-7° 7 7 5 commutes with the Dirac operator in (0) and decompose 



the Wigner function as 



18, 



■3 



*-ro s h =4(1 + /* (3) 



where we have omitted the flavour indices, k = k/|k| and a^ 1 , (p = 0,1,2,3) are the Pauli 
matrices, and h = ±1 are the eigenvalues of h. We multiply (J2J by pP, take the Dirac trace and 
integrate the hermitean part over ko- Introducing the 0th momenta of g„h, f^h = f(dko/2ir)guh, 

fij 



we note that the functions / i explicitly read 



ff h (x,k) = -j =! / ci 4 re"-''(^(x- r /2) 7 > K (x + r/2)), (4) 



dko 

~2~K J 


/" d 4 re ifc 


r ($hj(x 


-r/2)j°i> hi (x + r/2)), 


dko 
~2n J 


/" d 4 re ifc 


r ($hj(x 


-r/2)tl> hi (x + r/2)), 


dko 
~2n J 


/" d 4 re ifc 


r (4>hj(x 


-r/2)(-i 7 5 )^i(x + r/2)), 


dko 
~~2tt J 


/" d 4 re ifc 




-r/2) 7 7 5 ^(^ + ^/2)). 



Therefore, the /^(x, k) can be interpreted as follows: foh is the charge density, foh is the axial 
charge density, and fih and f 2 h correspond to the scalar and pseudoscalar density, respectively. 
From the Wigner-Dirac equation (J2J), we can now derive the following system of equations: 

f' Qh + i [M H , fih] + i [Ma, f2h] =0 , (5) 
f[ h + 2h\k\f 2h + i [M H , foh] ~ {M A , f 3h }=0 , 
f' 2h - 2h\\s\hh + {M H , hh} + i [M A , foh}=0 , 

f! 3h -{M H J 2h } + {M A ,fi h }=0, 

where the prime denotes a derivative with respect to r/. It is understood that M and the f^h 
are flavour matrices. Note that the commutators in l(HJ). which mix particle flavours, are essential 
for the production of the charges foh, and thus for our scenario. Moreover, one can infer, that a 
necessary condition for fL ^ is a nonsymmetric M. We already anticipated this when noting 
that for such a mass matrix the orthogonality of particle and antiparticle modes is violated. Note 
that the tree level dynamics given by Eqns. © closes for f^. When rescatterings, as described 
through nonlocal quantum loop corrections, are included, off-shell effects may become important, 
and one would have to solve for the full dynamics of the g^h- I n the nonrelativistic regime and 
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close to equilibrium however, in which off-shell effects are suppressed, it is possible to include 
rescatterings into Eqns. ® and still retain closure for the equations for the f^h- 

Now we fix the initial conditions for a Universe, which is void of fermions at the end of inflation. 
For an initially diagonal slowly evolving mass matrix, the Wigner functions for a zero particle state 



zero particles are (c/. Ref. 



foh - L h ' L h + R h ' R h' fih - - 2 ^( L h R h ab ), 

ff h = LfL\ - KRl ff h = 2Z(L* h a R b h ), (6) 



with 



T ab r Ug + hk b M* a 

L h — a b\ — 7T. ' K h ' 



1io a ' y/2u a (uj a + hk) 

where uj a = \/k 2 + |Af aa | 2 . Note, that for the case of real M aa , this reduces just to the usual 
choice of the components of the basis spinors in chiral representation. For a nondiagonal, but 
hermitean, M, one obtains the initial conditions by an appropriate unitary transformation. If 
additionally Ma / 0, as is the case for the SO(10) example discussed in the following, a biunitary 
transformation is necessary for diagonalization. 

Since foh is the zeroth component of the vector current, the charge of the species a carried 
by the mode with momentum k and helicity h is simply q a h{k) = f$h — 1- Note also, that the 
Lagrangean 

C = 4>aHa ~ MMH + il h M A )ba^a 

is f7(l) symmetric, and thus J2 a Qah(k) is conserved, as we shall verify explicitly for the SO(10) 
example discussed here. 

The scenario for coherent baryogenesis is as follows: initially, there are zero fermions described 
by appropriate initial conditions for the fah-i an d M is approximately constant in time. Then 
a phase transition occurs during which M changes rapidly, which leads to fermion production. 
Eventually, M stops evolving and the produced number of fermions as well as the charges /o^ 
stored within the different species are frozen in. We emphasise that /qI should not be confused 
with the number of produced particle pairs at preheating, which in our language can be expressed 
in terms of the (i = 1, 2, 3) as given in Ref. Q|. 
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3. THE BARR-RABY MODEL 

One possibility to break SO (10) down to the Standard Model is to use a Higgs multiplet A in 
the adjoint representation 45 and another pair of Higgses C and C in the spinor representations 16 
and 16. The apparently most simple implementation of this pattern of symmetry breaking, which 



is in accordance with particle physics observations, has been suggested by Barr and Raby |20j. In 
the following, we review the features of this model as far as they are relevant for our baryogenesis 
scenario, in particular we present in some detail the derivation of the higgsino-gaugino mass matrix. 
In the appendix, we give account of the conventions we use, in particular how the charges under 
the Standard Model group 

G SM = SU(3) C x SU(2) L x U(l)y (7) 

are assigned to the various multiplets of SO(10). 
We consider the superpotential 

WDkS(CC - /x 2 ) + ^- tr ^ 4 + ^M A trA 2 + T 1 AT 2 + M T T 2 2 (8) 

PA PA 

M s Ms 

where the additional fields S, P, Z\, Z 2 are singlets, T\ and T 2 10-plets of SO(10). Furthermore, 
there are the spinor C and the conjugate spinor C' . 

Let us first discuss the purely adjoint sector. The potential is at its minimum, when the 
condition 

dW 

-n = w = o (9) 

is met. When (A) = diag(ai, a 2 , a>3, a^, 05) &> i<72, it follows 

-^-af-M Aai = 0. (10) 



This can be solved by either ctj = 0, or = a, where 



MaMs nn 
a = ±y- ; - . (11) 

In order to step towards the Standard Model, it is possible to break SO (10) down to the left-right 
symmetric group 

G LR = SU(3) C x SU(2) L x SU(2) fi x U(1) b _l (12) 
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by choosing for (A) the Dimopoulos-Wilczek (DW) form 




a 



(A) 



a 



(13) 







V 







/ 



Note, that (A) being of DW form is proportional to the [B — L) operator given in Eqn. (1X261) . In 
the appendix, we give account of the explicit construction of the tensor- and spinor representations 
of SO (10) and the conventions we use. 

The two Higgs doublets of the MSSM are contained within T\ and are identified with the four 
components which remain massless by the superpotential (jHJ) when using the DW-form for (^4) . The 
additional six degrees of freedom of Ti, two colour triplets, become heavy and hence invisible at 
low energies, such that there is doublet-triplet splitting. The second 10-plet T 2 becomes necessary 
since a direct mass-term for the triplet components of T\ would lead to a disastrous rapid higgsino- 
mediated proton decay. 

The Higgs fields C and C reduce the SO(10) symmetry to SU(5). When minimizing the scalar 
potential, the absolute values of their VEVs are (C) = (C) = fA, and they point in the SU(5)-singlet 
direction with the quantum numbers of a right-handed neutrino. 

Both sectors, the spinorial and the adjoint, in combination reduce the SO(10)-symmetry to the 
Standard Model group Gsm- However, they need to be linked together in order to get a congruency 
of the assignment of Standard Model quantum numbers and to remove all pseudo-Goldstone modes 
from the particle spectrum. The obvious candidate term to add to the superpotential, CAC, 
however destabilizes the DW form (|13|) by altering the expression for the .F-term ()1U|) when the 
spinors get a nonzero VEV. Barr and Raby therefore suggested to add the additional spinors C 
and C' , which get a zero VEV. The conditions for potential minimization now become 



When comparing with Eqn. (|A26|) . we note that in the DW-form (|13[) we can identify (A) = 




(14) 



(15) 



§ia(i? — L). If we assume that the VEV of P is fixed, then Z\ and Z 2 settle to 




(16) 



(17) 
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since C and C point into the right-handed neutrino direction, where B — L = 1. We have hence 
achieved a link between the spinorial and adjoint sector without changing the form of — F\. 

In our model, CP-violation arises from the phase between £ and £ and therefore from couplings 
of the adjoint to the spinor multiplets. Let us label the multiplets of the Standard Model group (JJJ) 
by K. The representations 16 and 45 harbour as multiplets with common Gsm quantum numbers 
K = (3,2, g), K = (3,1,— |) and K = (1,1,1). These multiplets therefore mix through the 
higgsino mass matrix. The corresponding conjugate multiplets in 16 and 45 are labeled by K. 
Furthermore, all these representations contain the singlet (1,1,0). 

The spinor pair with 32 complex degrees of freedom breaks the 45-dimensional SO(10) down 
to the 24-dimensional SU(5). The 21 Goldstone modes come from the multiplets K = (3,2, g), 
K = (3, 1, — |), K = (1, 1, 1) plus one linear combination of the singlets K = (1, 1,0) within 16 
and 16. The 45-dimensional adjoint reduces the SO(10)-symmetry to the 15-dimensional Glr- 
Because of the DW VEV being proportional to the (B — L) operator, the 30 Goldstone modes can 
be identified with the multiplets for which B - L / 0, that are all colour triplets. 

Hence, by the supersymmetric Higgs-mechanism, there is a mixing of the higgsino modes with 
the gaugino sector, through the Lagrangean terms 

V2gv*T a i[)\ a + h.c. , (18) 

where T a is a generator of SO(10), normalized as tr(T a ) 2 = 1, A a a gaugino and <~p the scalar 
superpartner of the ^-fermion. This will induce higgsino-gaugino mixing mass terms for both 
multiplets, A and C when K = (3, 2, g) or K = (3, 1, — §), and only for C, when K = (1, 1, 1). 

Let us consider possible mass terms involving only the adjoint Higgs. If we denote the compo- 
nents of either (3, 2, g) or (3, 1, — |) by bx, we have (cf. appendix 

tiA 2 = -6a 2 -2b K b R , (19) 
tr^ 4 =6a 4 + Aa 2 b K b R + b 2 K b 2 R . (20) 

Hence, the portion ^p^-tr^ 4 + ^M^trA 2 of the superpotential (jSJ) gives for these modes a zero mass 
term 

™k = ^~M a = for if = (3,2 ,h and K= (3,l,-~), (21) 
Ms o 3 

where we have used the VEV (|11|) for a. This result is expected, since the multiplets in question 

are Goldstone. In contrast, we find 

CtCL 2 

m K = — for K = (1,1,1), (1,3,0), (22) 
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m K = -2— for if =(8,1,0). (23) 
M s 

Let us now discuss the mixing of the adjoints and spinors. ipA K and ipQt_ get a mixing mass 
term through 

SA R 5C' R ^2M S 

The derivation of this term is instructive and works as follows: In the block-diagonal basis, we 
have 



A W = I _ I ' (25) 



which transforms according to (jA30|) to the off-diagonal basis as 

1 




A — 7T- 1 A BL TT - 

^10 ~~ U BLOCK /1 lQ U BLOCK — ^ 



' 10 ilO 



ilO -10 



(26) 



The single degrees of freedom 10 are represented by the ten antisymmetric 5x5 matrices with two 
nonvanishing entries of the value 1 / y/2. Without loss of generality, we pick the matrix with — 1/ y/2 
in the first row, fourth column, corresponding to one degree of freedom of the K = (3,2,— g)- 
multiplet. In order to let the tensor A act on a spinor, we make use of Eqn. (jA24|) and represent 
it in terms of the T-operators as 

A = - ( [Tt , r 4 ] + i [r x , r 9 ] + i [r 6 , r 4 ] - [r 6 , r 9 ] ) . (27) 

lb 



When paired with the spinor 



1 r 



2 



X3X2X5 



10), (28) 



a Gsm singlet is formed, and after anticommuting the Xi-op era t° rs > we find 

(0|xiX2X3X4X5^ = ^ , (29) 

from which we immediately obtain Eqn. ()24|) . The higgsino-gaugino mixing terms can be derived 
from (|18f) in a very similar way. 

Finally, for the mixing among the spinors we have 

PW (30) 



5C R 5C' K s A M s ' 
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where olk = |i [(B — L)k — 1] or explicitly, 



OLK 



-1 for K = (3,2,1) 
-2 for K = (3,1,-f) 
for K = (1,1,1) 



(31) 



We have used here the VEVs (|16M17|) and again the proportionality of (A) to the (B — L) operator. 
We are now in the position to write down the higgsino-gaugino mass matrix: 



4>X K ipA K ipc K fa', 

( -iV2 7K ga g(C) 
iV2>jKga m K 




(32) 



V 



g(c) 
o 



o \ 

■ (C)(P) 
V2M S 
a(P) 



k(S) ia K ^-j^ 



( (C) (P) 
^ V2M S 



>«(P) 



M a J 



+ h.c. 



where 



IK 



\ for ^=(3,2,i) 
1 for K = (3,1,-|) 



(33) 



for K = (1,1,1) 

The mass matrix is nonsymmetric, therefore being endowed with the necessary prerequisites for 
coherent baryogenesis. 

4. SIMULATION OF COHERENT BARYOGENESIS 



The superpotential is of the type suitable for hybrid inflation. We assume that symmetry 
breaking by the adjoint sector has already taken place before or during inflation and is preserved 
throughout the subsequent history of the Universe, such that possible monopoles are diluted. We 
therefore do not consider the dynamics of the field A. For a discussion of the role of cosmic strings 



formed at the transition Glr — > G$m after inflation, we refer to Ref. 
The part of the scalar potential relevant for hybrid inflation reads 

2 



V = K A 



+ 2k z \SC\ 



(34) 



where we have used C = C* due to the vanishing of the D-terms and have written C = C. 

During inflation, the VEVs of C and C are sitting at a minimum located at zero, and S rolls 
down a logarithmic slope until reaching the critical value S cr = /J,, such that the value zero for C 
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and C becomes a maximum. The waterfall regime begins, at the end of which the scalar fields 
settle down to the supersymmetric (V = 0) minimum S = 0, \C\ = \C\ = [i. This is a rapid phase 
transition which brings coherent baryogenesis along. We simulate this scenario for the parameter 
k = 0.05, a damping rate T = 0.02fi and also take account of the expansion of the Universe, cf. 

1.5LA 



0.5LA 



-0.5LA 




10 20 30 40 50 60 70 80 90 100 

FIG. 1: Epoch of phase transition in the SO(10)-modcl 



Damping partly comes into play because of the perturbative decay of the inflaton. More impor- 
tant at the beginning of the waterfall regime is however the phenomenon of tachyonic preheating: 
Since the fields C and C attain a negative mass square term, modes with momenta less than this 



lis effect b y in troducing the damping rate 



23 



281 ] . Damping will also 



mass get produced exponentially fast. Here we mimic t 
r. For numerical studies of this process, see Refs. | 
receive a contribution from fermionic preheating. A proper treatment of fermionic preheating, 
which includes rescatterings, would require techniques used in Ref. [3], which have so far not been 
applied to the question of inflaton thermalisation through decay into fermions. 

In order to keep the discussion simple, we do not take the dynamics of the singlet fields Z\ 
and Z2 into account here. In principle, their VEVs only get fixed when C and C acquire nonzero 
VEVs. A possible way to fix Z\ and Z2 already during inflation is for example to shift the spinors 
away from the zero VEV, as proposed in Ref. [5J and is also applicable to SO(10)-models p. 

For the remaining parameters, we choose fi = 0.5 x 10 16 GeV, Ms = 550/i, Mc = 0.02/i, g = 0.2, 
C = -0.02, f = 0.05i, a = 25/i and (P) = 24/Lt. 
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FIG. 2: The produced charges for the multiplet (3, 2, g). 
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FIG. 3: The produced charges for the multiplet (3, 1, — |). 



The charge numbers which are plotted in figures [21 an d El refer to the mass matrix, which is 
diagonal after the phase transition is completed. Hence, there is mixing among ip\ K , tpA K , ^c K 
and i^c'x- The mass matrix is diagonalized via a biunitary transformation. It turns out that 
the higgsino-gaugino mass matrix in the supersymmetric vacuum has two heavy and two light 
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eigenvalues. We only display the charge production corresponding to the light mass eigenvalues 
which we label by q\ and c/2- There is no substantial charge asymmetry stored within the heavy 

4 

flavours, |<?3|, |(/4| <C |<Zi|,|<72|- The apparent symmetry Qi = results from the overall U(l) 

i=i 

symmetry of the fermionic fields and is a useful check for the numerical results. We shall assume 
here that the resonant decay into gauge bosons and scalar Higgs particles is not important. This 
can be justified by noting that the necessary conditions for a resonant decay into bosons are not 
met for our choice of parameters. 

The qi are charges stored within the mass eigenstates which are Dirac fermions of the generic 
mixing form 

/ 



i a \ K X K + ia AK i>A K + ia CK ^c K + i^c'J>C K 



1,2,3,4, (35) 



where the coefficients ia\ and ij3yr are determined by the biunitary transformation diagonalizing 
the higgsino-gaugino mass-matrix. For the example we discuss here, the transformation matrices 
and the coefficients are determined numerically. 

Now, we discuss how the qi charges get transformed to B — L charge stored within fermionic 
matter. Decay into SUSY Standard Model particles can take place through the gaugino components 
Xk, A/f and through the Higgsinos ipc R > rfc-' ^he relevant operator for gaugino decay from the 
gauge supermultiplet Lagrangean is 

V2g {XFF + h.c.} (36) 

while the higgsinos decay through the dimension five couplings 

CFCF 



< 37 > 

CT a FCT a F . . 

"» Ms ■ (38) 
■73 CT ° r ^ rartf , (39) 

which are added to the superpotential (jSJ) , and where F are the Standard Model matter fields and 
the right-handed neutrino, contained in 16, and T denotes the operators defined in (|A6|) and (|A7|I . 
The index i = 1,2,3 denotes the matter generation. 

While the coupling ()36|) is universal for all three generations of Standard Model matter, 
*72 and *73 may be different for the three generations. We assume that 3 72.3 3> 2 72,3 S> 1 72,3 ; such 
that only the 3 72 3 are of relevance for the decays. In contrast, we require that the decays through 
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2 7i and 3 7i are not possible since the corresponding right-handed neutrinos are heavier than the 
decaying particle, such that only 1 ji is relevant. We now argue of which order these couplings 
should be for realistic scenarios. 



For nonthermal leptogenesis 



33| . one usually assumes that two of the three Majorana neutrinos 



from the different generations of matter fermions are heavier than half of the mass 

mj = y^K/i (40) 

of the inflaton fields, which are the ^ c -like components of C and C and the singlet S, such that 
their decay into right handed neutrinos is kinematically forbidden. Through the coupling Q37JI the 
right handed neutrionos acquire Majorana masses 

such that the requirement 2 ' 3 m u c > mi/2 reads 

where we have used (C) = \i. It appears to be reasonable that also the couplings 2 ' 3 72,3 are of the 
same order, as we shall assume. 

For the scenario we discuss, the lightest right-handed neutrino is important for the reheating 
process. The coupling (|37[) also allows for the decay of the inflaton fields at the rate 

The Universe becomes radiation dominated and entropy production stops, when T v = H, where 
H denotes the Hubble expansion rate. The reheat temperature at this time is 



T R = 0.555* Vr>npi , (44) 

and we take the estimate = 220, the number of relativistic degrees of freedom after reheating. 

The value for the Planck mass is mpi = 1.2 x 10 19 GeV. The mass of the lightest right-handed 

neutrino is therefore proportional to the reheat temperature: 

i 

W = 7.7x4=— 2«. (45) 

Taking for the highest reheat temperature allowed by the gravitino bound Tr = 10 11 GeV and 
for the parameters we use, we find 1 m u c < 5 x 10 9 GeV while 2,s m u c > 2 x 10 14 GeV. Therefore, 



15 



a fortuitous hierarchy of five orders of magnitude for the Majorana masses is required. This is 
usually assumed for scenarios of nonthermal leptogenesis. 

For the coherent baryogenesis model we consider here however, we can also allow for all Majorana 
masses to be larger than mj/2. Under these circumstances, S and the ^ c -like components of C and 
C cannot decay through the term (|H7jl into two right-handed neutrinos but decay to three particles 
via the operators (|38|) and (|39|) . Since these processes involve no Majorana particles, leptogenesis 
is absent for this scenario. 

We also have to deal with the fact, that for ax = 0, namely for K = (1, 1, 1), ipc K an d ^C K do 
not mix, cf. the mass matrix (|32jl . Therefore we assume that also the fields C and C' may decay 
through couplings of the above type, suppressed however by additional powers of (R)/Ms, where 
R is some singlet with a VEV. 

Note however, that for K = (1,1,1) the mass matrix (|32j) is block-diagonal, such that only the 
pairs \k-iPc k an d ipA^-^c' are mixed. Only for the second pair, the CP-violating parameters £ 
and C are relevant and an asymmetry is generated, which vanishes however after the decay into 
matter. 

The charges qi gets processed differently to B—L when Vl/j decays through its various components. 
Let us denote the B — L number resulting from the decay of a component X of a *$>i quantum by 
Tx- By the couplings (jHEJi the reactions 

A*-— i^ + i/* (46) 
*k^Fk + " c * (47) 

are induced, where one of the particles on the right hand side is a scalar, the other one a fermion. 
Due to its Majorana mass term coming from the operator ()37|) . the right handed neutrino v c is its 
own antiparticle and therefore carries effectively B—L = at tree-level. The resulting B— L-charge 
is therefore the one stored within Fj> and F% and we find 

T Xk = 1> K=(3,2,±), (48) 



T Xk =-~, if =(3,1,--), (49) 



1 * = (3,1,-| 

T Xk =1, K= (1,1,1), (50) 

T ^k = -\> # = (3,2,—), (51) 

7a a =^, # = (3,1,1), (52) 

T Xr = -1, K = (1,1,-1). (53) 
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Similarly, the coupling (|3*7)) allows for the decay reaction 

fcfe— ** + "°*. (54) 

Hence, the charges get transformed to 

T c k = -\i K= (3,2,-±), (55) 

T c R = \. K=(3,l,^), (56) 

T Ck = -l, K = (1,1,-1). (57) 

We can calculate the term (|HH|) oc 7? using the techniques explained in the appendix EI It is 
however easier to note that (c/. Ref. [3oJ]) 

(3,2,i)®(3,l,i) D (l,2,i) C 10, (58) 

as well as 

(1,1,0)® (1,2,-i) = (1,2,-^) C 10. (59) 

The components of ipc K f° r ^ = (3, 2, g) therefore decay to 

i>c d ^d c * + e*, (60) 

^^rf c * + ^*, (61) 

where i/jq denotes the d-quark like higgsino, ipc u the it-quark like one. The charges hence get 
transformed as 

T Ck = 1, for IT = (3, 2, i). (62) 

The u c -quark like higgsino with K = (3, 1, — |) decays through the 73-coupling (|39|) . We note 
(cf. Ref. |3o|) 

(3, 1, -|) (3, 1, i) = (3, 1, -i) (6, 1,-1) C 120 , (63) 

(1, 1,0) ® (3, 1, i) = (3, 1, i) C 120 , (64) 

and therefore have the reaction 

i>c u c — d c * + <P , (65) 

and the charge conversion 

T Ck = \, for K = (3,1,-|). (66) 
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Finally, the e c like higgsino K = (1, 1, 1) turns into matter via the 72-coupling (|38[). as can be 
seen by 

(l } l } l)®(l,2,-~) = (l,2,i)cl0, (67) 
(1,1,0)® (1,2, -i) = (1,2, -i) C 10. (68) 



Consequently, the decay reaction is 



^C' ~^e* + v\ (69) 



and the resulting asymmetry 



T Ck = 2, for if = (1,1,1). (70) 



The procedure to obtain the produced B— L-density is as follows: We first integrate the produced 
charge numbers qi over momentum space in order to obtain charge densities Q{. From the biunitary 
diagonalization of M in the supersymmetric minimum the contributions of Xk, Xg, iftc K an d ipg- t° 
the *S?i are determined, which gives the branching ratios and therefore the respective contributions 
for the decays of the Dirac fermions to Standard Model matter. As a formula, this reads 



,0 



11 k 

i=l 



(71) 



where the factors of three come from the presence of three generations of matter and 



2 for K = (3,2,1) 

3 for K = (3, 1,-|) • (72) 
2 for if =(1,1,1) 

The total (B — L) number density produced at the phase transition is taking account of the 
multiplicity of colour and flavour given by 

n° B -L = 6n (3,2,I) + 3n (3,l,-|) • (73) 

A study of the parametric dependence of the produced asymmetry is beyond the scope of this paper, 
which is to show that coherent baryogenesis is viable with the gauge group SO(10), and shall be 
discussed elsewhwere. Therefore, we content ourselves with presenting just one typical numerical 
example here. The parameters yet to be specified are the 3 72,3 and 71. We can set effectively 
*7i = because it is either very small due to the restrictions given by the reheat temperature (|44|) 
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and the gravitino bound or, in the case of 1 m u c > mj/2 decays through the coupling Q37JI do not 
take place. In accordance with the relation (|42|) we choose 3 72 = 3 71 = O.O5M5//X. Then, we find 

n^_ L = 2.5 x l(rV . (74) 

In order to estimate the baryon to entropy ratio, we express the entropy density s through the 
reheat temperature Tr as 

s = 27rVTl/45, (75) 

and the Hubble expansion rate is given by 

H = 1.66^57— , (76) 
mpi 

where mpi = 1.22 x 10 19 GeV is the Planck mass. 

During the epoch of coherent oscillations, that is between the end of inflation and the onset of 
radiation era, the Universe is matter dominated and expands by a factor 



2 

a ( Ho \ 3 



where Hq is the expansion rate at the end of inflation, given by 



(77) 



Ho = J 8 i— 2 - (78) 
w 3 mp\ z 



Putting everything together, we find 

s 3s V a J 4 Vo 
where we have taken account of a division by three for sphaleron transitions promoting {B—L) to 
B asymmetry. 

The value for the vacuum energy at the end of inflation is Vq = k 2 /j, 4 , and by Eqn. (|79j) . we find 

^ = 1.0 x 10- 10 , (80) 

where we have chosen Tr = 2 x 10 10 GeV. Hence, it appears that in order to get a BAU in 
accordance with observation, there has to be a reheat temperature of order of the upper bound 
allowed by the requirement that gravitinos shall not be overproduced. However, our estimate of 
entropy production is rather crude. It is conceivable that initially the decay of C, C and S is 
enhanced due to tachyonic and also fermionic preheating. This could lead to an initial radiation- 
like equation of state (23] or shorten the matter-dominated era and would therefore lead to less 
dilution of the initial asymmetry. 
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For the case 1 m Uc < mi/ 2, it is of interest to compare the result (|8U|) with the baryon asymmetry 
resulting from nonthermal leptogenesis, which is given by Q 

^=0.5^, (81) 
s mj 

where 

is the maximal CP- violation which may arise from the decay of the lightest right-handed neu- 
trino |3J|, and m Ua denotes the heaviest mass eigenvalue of the light neutrino mass matrix. When 
we assume m V3 = 50meV and use the same parameters as for the coherent baryogenesis example 
and Eqn. ()45|) for 1 m uc , we find tib/s = 6 x 1CU 12 . Therefore our example corresponds to a point 
in parameter space where coherent baryogenesis dominates over leptogenesis. However, we expect 
that also the opposite case may occur for a different set of parameters. 

5. CONCLUSIONS 

In this paper, we show that during the phase transition terminating SUSY SO(10) hybrid 
inflation, a charge asymmetry within the Higgsino sector may be produced through the mechanism 
of coherent baryogenesis and subsequently turned into baryons. CP- violation is provided by the 
couplings of the spinorial to the adjoint representations and occurs at tree-level. This is very 
different from leptogenesis, a one loop effect, where CP-violation is sourced by the matrix of 
Yukawa couplings between the neutrinos and the Standard Model Higgs field. Since the spinor- 
adjoint couplings are an indispensable part of the Barr-Raby model, coherent baryogenesis naturally 
occurs at the end of SUSY-SO(IO) hybrid inflation. Together with similar results which we found 
for the Pati-Salam group [7j, this indicates that effects from fermionic preheating are generically 
of importance for the generation of the BAU in hybrid-inflationary scenarios. 

Coherent baryogenesis however has been neglected in the standard picture of baryogenesis in 



SUSY-GUT hybrid inflation so far, which is as follows |33j: The inflaton decays into right-handed 
neutrinos which then decay out-of equilibrium into Standard Model matter, leaving behind a 
B—L asymmetry via the leptogenesis mechanism 34]. Since the decaying Majorana neutrinos are 
not produced by the thermal background, this scenario is often called nonthermal leptogenesis. 
However, it imposes strong constraints on the hierarchy of the masses of right-handed neutrinos, 
as we discuss in section 0] We also emphasize that coherent baryogenesis relaxes this constraint 
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and allows for all Majorana masses to be larger than the inflaton mass, since the mechanism 



does not rely on leptogenesis and the decay of Majorana particles. In conclusion, the re 



ations of 



the parameters of hybrid inflationary models to the BAT 1 as suggested c.e. in Rols. jj], Q] by 
considering leptogenesis, should be altered, since coherent baryogensis turns out to be an additional 
source of the BAU, which may dominate in some regions of parameter space. 



We emphasize that nonthermal leptogenesis and coherent 



34 



jar yog enesis should not be confused 
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36 
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3a ] . an appealing feature 



with the often discussed thermal leptogenesis mechanism 
of which is that the BAU is generated from a Universe which is - within horizon scale - in thermal 
equilibrium. Leaving aside primordial density fluctuations, all cosmological observables including 
the BAU would then be predictable from an effective theory valid up to the Majorana mass scale 
of the neutrinos. 

Grand Unified Theories open up many possible paths for the generation of the BAU and it is yet 
not known where the actual asymmetry originates from. The various mechanisms allow to establish 
relations to the paramters of the underlying models. While leptogenesis renders constraints on the 
neutrino sector, coherent baryogenesis is of interest since it is a scenario of GUT-baryogenesis and 
thereby related to the dynamics of symmetry breaking. 
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APPENDIX A: SO(10) 

Besides by tensors, orthogonal groups may also be represented by spinors, which satisfy a Clif- 
ford algebra. In order to construct group-transformation invariants, both types of representations 
need to be linked together via Dirac gamma matrices. For the familiar case of the Poincare group 
SO(3, 1), it is often convenient to use a specific representation for these matrices. In contrast, 
one better circumvents the tedious task of explicitly constructing ten 32 x 32 gamma matrices for 



SO (10). Mohapatra and Sakita |39| have therefore devised a very useful technique for performing 
calculations involving spinors and tensors, employing just abstract commutation and anticommu- 
tation relations. 

On the other hand, when it comes to symmetry breaking, one has to choose a certain convention, 
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that is a certain basis, how the particles of the Standard Model are assigned to the representation 
16 of SO(10). This assignment fixes in turn the definition of the charge operators and hence the 
quantum numbers of certain entries in vectors and tensors of SO(10). 



While the paper by Mohapatra and Sakita 



39] does not provide much details of tensor rep- 



resentations and symmetry breaking, such a discussion can be found in the comprehensive work 
by Fukuyama et al. jsQ], where in turn spinors are neglected. The coupling of spinors to ten- 
sors is explained for SO(10) by Nath and Syed j^. In the paper by Barr and Raby [^(J, which 
contains the model we consider here, a basis where tensors nicely decompose into blocks of SU(5)- 
representations is chosen. Unfortunately, the choice of basis and normalizations is not explicitly 
given, but has to be inferred by the reader. 

In the following, we give some detailed account of the construction of SO(10)-singlets, follow- 
ing the conventions of Barr and Raby. Explicit expressions for the charge operators acting on 
spinors and tensors as well as for the accommodation of the Standard Model particles and the 
right-handed neutrino in the representation 16 are given, which shall ensure an easier and faster 
comprehensibility of the Barr and Raby analysis as well as of our calculations. 

Charge Assignments 

We denote by Q the electric charge, by Y the weak hypercharge and by l\ the weak isospin. The 
charges which are not gauge symmetries of the Standard Model are baryon minus lepton number 
B — L as well as the SU(2)^-isospin 1^ and the less known charge X. There are linear dependencies 
among these charges, which are given by 



Q=ll + Y, (Al) 

B - L=2(Y - I 3 R ) , 

B - L=-UY - X) . 
5 

Note that, when comparing to the conventions by Fukuyama et. al. |3o| . we have twice as large 
values for (B — L), such that for a single lepton, we have (B — L) = — 1. 

In table E we give the charge numbers of the Standard Model particles and of the right-handed 
neutrino. 
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TABLE I: Quantum numbers of matter 












Y 


B - L 


X 


Q 


M 


2/3 


1/2 





1/6 


1/3 


-1 


_ U 








\ / 


-1/3 


-1/2 





1/6 


1/3 


-1 




u c 


-2/3 





-1/2 


-2/3 


-1/3 


-1 




d c 


1/3 





1/2 


1/3 


-1/3 


3 


L 







1/2 





-1/2 


-1 


3 


■(:) 














-1 


-1/2 





-1/2 


-1 


3 












-1/2 





1 


-5 




e c 


1 





1/2 


1 


1 


-1 



S0(2A) in an SU(A) Basis 



This section contains a review of the paper by Mohapatra and Sakita 39], but adopts the basis 
conventions of Barr and Raby [^(J. 

Let us introduce iV operators Xi (* = l>---)-^0> acting on an antisymmetric Fock space, which 
obey the following anticommutation relations: 

{^,x]}=Sij, (A2) 
{Xi,Xj}=0. (A3) 



The operators defined as 



satisfy the SU(A r ) algebra: 



T) = xlXj (A4) 



[T i j ,T k l )=5 k j T l l -5\T k j . (A5) 

We now introduce the 2A^ operators 

^3=-KXj - x}) , j = l,..-,N, (A6) 

T N+j = Xj + x], (A7) 

which obey by Eqns. (jA2| IA3|) the Clifford algebra 

{T i ,T j } = 25 ij , i,j = l,...,2N, (A8) 
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and hence, the algebra of generators of S0(2iV) is given by 

^ij = ^i,^\- (A9) 

Since the dimension of the spinor representation of SO(2iV) is 2^, a concrete representation could 
be constructed for SO (10) in terms of 32 x 32-matrices, which however shall not be done here. 

The spinor states can be constructed by letting the N creation operators x\ ac * on the "vacuum" 
|0), such that the spinor representation is 2^-dimensional, as it should. 

It is well known, that the spinor representation of SO(2iV) is reducible. We therefore define 

2N N 

r = i 7V n r « = IK 1 - 2n i)' ( A1 °) 

i=l j=l 

where we have introduced the number operators 

nj = X)X3 ■ (AH) 

The chiral projectors ^(1 db Tq) give therefore rise to the two irreducible 2^" ^dimensional repre- 
sentations containing only even (case "+") or only odd (case "— ") numbers of creation operators. 
Now let ^ be an SO(2N) spinor state. We are interested in calculating products of the form 

* T BT h ...r iM y, (A12) 

involving a certain number of T matrices. The matrix B is necessary since ^ T does not transform 
as a conjugate spinor when acted upon with an infinitesimal SO(10)-transformation ey: 

(J^ieyEij-*, (A13) 
^t=-ie ij #E ij , 



We require from B the property 



5* T =i* T eyE£. . 



B-^B = -Ey , (A14) 



such that 

5(V T B) = ie^BB-^B = -ie i3 {^ 7 'B)Ey , (A15) 
i.e. ^ T B transforms as a conjugate spinor. The condition (|A14|) can be met if 

B^YjB = ±Ti . (A16) 
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By choosing the minus-sign in the latter equation, we find 

N 

£=n r *> ( ai7 ) 

i=i 

because for i = 1, N the Tj are represented by antisymmetric matrices, while for i = N+l, 2N 
by symmetric ones. 

For N = 5, we can arrange the Standard Model particles in the spin-16 representation, which 
is projected out of the 32-dimensional spinor * by i(l-r )*. Defining 

u* = \e lkm x\x\xl\V) , (A18) 

^=^ H45 xlxMlo>, 

u C i=xUlxl\0) , 
dt=x\\0), 

f=xJ|o>, 



e=xS|0>» 

^ c =x{x2xtxlxt\o) 

e c =x\xlxt\0) , 



where i,k,l = 1, 2, 3. Cf. also Ref |41|], where the doublet and triplet blocks are interchanged. 

The next task is to construct the charge operators. For example, the ladder operators associated 
with the left isospin take u d and v <-» e. They are therefore given by 

It=x\x±, (A19) 

^Z=xlx5 • 

The weak isospin operator is hence 

lL = ~[liJZ} = l(n 5 -n4). (A20) 

By comparison with the charge numbers in table [IJ we can identify 

3 5 

Y = 3 - \1L n i = 12- d r i' r 6] + I r 2>r 7 ] + [r 3 ,r 8 ]) - - ([r 4 ,r 9 ] + [r 5 ,r 10 ]) , (A21) 

i=l j=4 

where we have used 

[T 5+j ,T j ] = -4in j + 2i. (A22) 
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When identifying the indices of the T operators with matrix rows and columns as implied by 
Eqn. (jMJ), we can explicitly write down the suitably normalized Y in tensor representation: 

Y = diag (1/3, 1/3, 1/3, -1/2, -1/2) ® a 2 . (A23) 



Generally, we use as rule for conversion of the operator to the tensor representation 1 

"1 

_4 

which reads for the special case of the charge operators 



= $ia,5jb ~ Sj a 6 ib , (A24) 

ab 



^ [r 5+i ,ri] = m- - = diag^ij,^,^,^,^) ® o 2 . (A25) 

Now, we easily find the other charge operators. Putting everything together, we have in spinor 
and in tensor representation 

1 3 

g=-^n«-n4 = diag (1/3, 1/3, 1/3, -1, 0) ® a 2 , (A26) 
6 i=i 

ll=\inh-n 4 ) = diag (0,0, 0,-1/2, 1/2) ®a 2l 

4=^(1 - n 4 - n 5 ) = diag (0,0, 0,-1/2, -1/2) ® (T 2 , 

2 3 

B-L=-Y i n i -l = diag (2/3, 2/3, 2/3, 0, 0) ® a 2 , 

1 = 1 

3 5 

Y =\ E n i - J £ n J = dia § (V3, 1/3, 1/3, -1/2, -1/2) 8) a 2 , 
i=l 3=4 
5 

X=-2 n i + 5 = dia g (- 2 ' ~ 2 > ~ 2 > ~ 2 > ~ 2 ) ® ^2 , 
i=l 

where we have used the normalization convention (|A1|) . 

The operator representation for the charge operators Q is suitable for finding the charge eigen- 
values q of the spinors through Q^f = q^f. 

Tensors can be constructed from the fundamental 10-dimensional vector $io by taking antisym- 
metric products, such that a rank n tensor is of dimension 10 • 9 • ... • (10 — n + l)/n\. Explicitly, 
for the vector and the rank two tensor, the charges implied by the gauge-covariant derivatives are 
given by the eigenvalue equations 

e,*io=g*io, (A27) 

[Qq, $45] =9$45, 

where Q is acting here by matrix multiplication. 



1 This is of course strictly speaking no equality but an assignment of an operator acting in Fock space to an operator 
acting in tensor space. 
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The Tensor Representations 

In order to perform calculations such as 16.45.16, tr45 4 and 16.10.16, we need to identify the 
Standard Model multiplets within 10 and 45, just as we did for the 16 in (lA"T8l) . We first note, 
that under SU(5), the fundamental representation of SO(10) decomposes as 10 = 5 © 5. Let us 
denote an element of 5 in the representation 10 of SO(10) by ^fo, an element of 5 by <3?f . Since 
they obey 



X$l = 2$f and X$l 



-2$ 



10 > 



(A28) 



they are of the form 



10 



1 

72 



a 5 
-iai 



and $fo = ^ 



a 5 
iai 



(A29) 



V - ia 5 / \ ia5 / 

with J2i=i \ a i\ 2 = 1 an d Ya=i \®i\ 2 = 1- To remove this inconvenient mixing of the upper and lower 



five-blocks, we introduce the unitary transformation 



U BLOCK - ^= 



1, il. 



BLOCK 




(A30) 



such that 



t^BLOCK^lO 



a 5 




/ n \ 



and U-i 



BLOCK^lO 



(A31) 



V / V «5 j 

By this change of basis, the charge operators become diagonal, for example 

/ 

^BLOCK^t^BLOCK = 2 




(A32) 
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We can therefore immediately see how the entries of 45 transform under SU(5), namely 

/ 

^BLOCK^S^BLOCK = 



24 1 


10 


10 


24 ©1 



(A33) 



where the single entries represent 5 x 5-blocks and the blocks in the upper left and the lower right 
are to be related to each other by the factor of minus one. The SU(5)-singlet 1 has here the form of 
the matrix l/\/5l5- The arrangement of the G^M-multiplets contained in 24 can be schematically 
written as 





12 3 


4 5 


1 

2 
3 


(8,1,0)0(1,1,0) 


(3,2,-|) 


4 
5 


(3,2,|) 


(1,3,0) 



(A34) 



and finally 10 of SU(5) decomposes into 





1 2 3 


4 5 


1 

2 
3 


(3,1,-|) 


(3 5 2,|) 


4 
5 


"(3,2, |) 


(1,1,1) 



(A35) 



where the matrix is imposed to be antisymmetric, since it is identified with the antisymmetric part 
of 50 5 of SU(5). 
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